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ABSTRACT 

First, we describe a general procedure to produce high quahty vector mag- 
netograms using the Imaging Vector Magnetograph (IVM) at Mees Solar Obser- 
vatory. Two IVM effects are newly discussed and taken into account: (1) the 
central wavelength of the Fabry-Perot is found to drift with time as a result of 
undiagnosed thermal or mechanical instabilities in the instrument; (2) the Stokes 
l^-sign convention built into the IVM is found to be opposite to the conventional 
definition used in the study of radiative transfer of polarized radiation. At the 
spatial resolution 2" x 2", the Stokes Q, U, V uncertainty reaches ~ 1 x 10~^ 
- 5 X 10~^ in time-averaged data over 1-hour in the quiet Sun. When vector 
magnetic fields are inferred from the time-averaged Stokes spectral images of Fel 
6302.5 A, the resulting uncertainties are on the order of 10 G for the longitudinal 
fields (-By), 40 G for the transverse field strength and ~ 9° for the magnetic 
azimuth (0). The magnetic field inversion used in this work is the "Triplet" code, 
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which was developed and implemented in the IVM software package by the late 
Barry J. LaBonte. The inversion code is described in detail in the Appendix. 

Second, we solve for the absolute value of the vertical electric current density, 
\Jz\-, accounting for the above IVM problems, for two different active regions. 
One is a single sunspot region (NOAA 10001 observed on 20 June 2002) while 
the other is a more complex, quadrupolar region (NOAA10030 observed on 15 
July 2002). We use a calculation that does not require disambiguation of 180° 
in the transverse field directions. The \ Jz\ uncertainty is on the order of ~ 7.0 
mA m~^. The vertical current density increases with increasing vertical magnetic 
field. The rate of increase is about 1 — 2 times as large in the quadrupolar NOAA 
10030 region as in the simple NOAA 10001, and it is more spatially variable over 
NOAA 10030 than over NOAA 10001. 

Subject headings: Sun: magnetic fields - methods: data analysis 



1. Introduction 



The Imaging Vector Magnetograph (IVM) was built in the early 1990s at the Mees 
Solar Observatory, University of Hawaii, to measure vector magnetic fields on the Sun. 
The IVM is a complicated instrument requiring care in both taking of the data and the 
subsequent calibration and analysis for the determination of accurate vector magnetic fields. 
Characteristics of the instrument and the data calibration have been described in a series 



of published papers (Mickey et al. 1996; LaBonte, Mickey &; Leka 1999; Leka & Rangarajan 



2001 LaBonte 2004). The present paper is a continuation of this series in which we identify 



and discuss two important issues with the IVM data and analysis that were previously not 
acknowledged in the literature. 

The IVM is designed to combine the advantages of two distinct types of magnetograph, 
namely, the narrow-band filter polarimeters and the spectroscopic polarimeters. Both types 
are based on the Zeeman effect with the use of magnetic-sensitive spectral lines. They both 



are in use, for example, onboard the newly launched HINODE spacecraft (Kosugi et al. 



2007). 



Narrow-band filter magnetographs employ birefringent filters for wavelength selection 



(Hagyard et al. 1982 Ai & Hu 1986 Zirin 1995). By rotating the crystal elements inside 



the filter, a narrow band spectral image of a polarization state is formed at a magnetically 
sensitive spectral line position. The advantage of the filter devices is that high spatial 
and temporal resolution 2-dimensional polarization images can be formed at high signal- 
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to-noise ratios. It takes only ~10 minutes to complete a single vector magnetogram. The 
disadvantage is that the four Stokes parameters are not determined simultaneously within 
this period, so that spurious polarization may be introduced due to changes on the Sun and 
in the Earth's atmosphere. Neither is the entire spectral line observed, which limits the 
diagnostic methods used in the magnetic field inversion process. The magnetic field filling 



factors cannot be estimated with a narrow-band filter magnetograph (Lites, Martmez Pillet, 



fc Skumanich|1994 ) and magneto-optical effects may be underestimated in the measurements 
of transverse fields (Hagyard et al. 2000). 



Spectro-polarimeters record the solar images using a narrow slit placed in front of a 
spectrograph. Stokes parameters are measured as functions of wavelength across a spectral 
line sensitive to the Zeeman effect. High spectral resolution can be achieved. To minimize 
the influence of time-dependent seeing and solar variations, one such spectro-polarimeter, 
the Advanced Stokes Polarimeter (ASP), uses two CCD cameras to record orthogonal polar- 



ization states simultaneously (Lites, Martmez Pillet, & Skumanich 1994 Lites 1996). The 



entire Stokes spectral profiles are taken into account for inferring the magnetic fields and the 
magneto-optical effect and magnetic filling factor can both be well estimated. A disadvan- 
tage is that this kind of instrument takes a longer time, perhaps ~20 minutes to an hour, to 



obtain a single vector magnetogram at high signal-to-noise ratios (Elmore et al. 1992 Lites, 



Martmez Pillet, & Skumanich 1994 Lites 1996). Furthermore, the observing field of view 



is often smaller than those of narrow band polarimeters. The Haleakala Stokes Polarimeter 



(HSP) also belongs to this category (Mickey 1985; Canfield et al. 1993) 



The Mees IVM measures the vector magnetic field in a slightly different way from the 
above instruments. A Fabry-Perot interferometer is employed to scan the entire magnetically 
sensitive spectral line (Fel A6302.5) at each imaging pixel with a spectral resolution of 0.04 
A and a wavelength coverage of 1.2 A. In principle, the IVM is in the spectro-polarimeter 
category but it has a wide field of view, 4.7' x 4.7'. The design aims at tackling the atmo- 
spheric seeing variation, at the same time acquiring data at high temporal (~ 2 minute data 
cadence) and spatial (0.55" x 0.55" pixel size) resolution. The polarization precision is aimed 



at 0.001 (Mickey et al. 1996) 



As pointed out by LaBonte, Mickey & Leka (1999), two basic steps are involved in the 



reduction of the IVM data. Step 1: deduce the Stokes parameters J, Q, U, V from the raw 
IVM data. The Stokes I represents the total intensity of the beam; the Stokes Q, U represent 
states of linear polarization in which the electromagnetic waves are in planes separated by 
45° with planes are perpendicular to the direction of the beam; the Stokes V represents the 
circular polarization state in which the electromagnetic waves rotate around and propagate 
along the direction of the beam. The IVM software package was written to deal with many 



-4- 



issues in the data in this step (Mickey et al. 1996 LaBonte, Mickey & Leka 1999 Leka 



& Rangarajan 2001; LaBonte 2004). Step 2: infer the vector magnetic field from Stokes 



parameters via the radiative transfer equations for polarized spectral light in the presence 
of magnetic fields. This step is dependent on Stokes inversion methods. Three methods are 
implemented in the IVM package: "wavelet method" , "derivative method" and the "triplet" 



method. A number of investigations based on IVM data used the first two methods [e.g.] (Leka 
fc Skumaniclil[T999| [Regnier et aT][2002| [Bleybel et aL][2002| |Leka fc Barnesl|2003a||bt [Welsch 



et al. 2004; Metcalf et al. 2005; Barnes et al. 2005 Leka & Barnes 2007). The current work 



will use the "triplet" method which will also be described in detail in the Appendix A. The 



"triplet" method was also used in the studies based on IVM data by a few [e.g.](Georgoulis 



& LaBonte 2004; Regnier & Canfield 2006 Georgoulis & LaBonte 2006 2007) 



Recently, two problems in the IVM data caught our attention. Neither of them has 
been discussed in the literature, nor have they been taken into account in the standard IVM 
package. The effects of these two problems are examined in this paper. The first problem is 
that the spectral line center produced by the Fabry-Perot shifts with time in the raw data. 
The cause of this problem is unknown, but is presumably related to a mechanical/thermal 
shift in the Fabry-Perot causing the spectral line to drift throughout the day. If uncorrected, 
the wavelength shift in 1 hr reaches about 0.2 A, equivalent to the (very significant) Doppler 
shift of 9.5 km s~^. We will discuss a procedure for minimizing the error caused by the shift 
in the Stokes data without jeopardizing the solar velocity measurements. We do not know 
for how long the spectral line drift has been a problem with IVM: conceivably it affects all 
data taken with this instrument. 

The second problem is that the sign convention built into the IVM determination of 
Stokes V is unfortunately opposite to the conventional definition employed in the analysis of 
the radiative transfer of polarized radiation. A brief check shows that the problem was not 
present in the data taken on 19 June 1998, but was present in data taken on 15 July 2002, 



and on 6 August 2003 (see Figure 5 of LaBonte (2004)). According to the (anonymous) 



referee of this paper, the problem has been unofficially discussed among IVM users, notably 
the late Tom Metcalf, who mentioned the Stokes \^-sign problem. By all the evidence, we 
suggest that the Stokes- \^ signs were correct before January 1999 with image size 256 x 256. 
The signs are wrong for the data taken after January 1999 when the IVM underwent a major 
upgrade, and acquired a larger CCD camera allowing image size 512 x 512. The Stokes V^-sign 
problem results in over-correction of the magneto-optical effect in the Stokes polarization 
measurements and must be removed if accurate inversions are to be obtained. 



An important application of vector magnetograph is to determine the vertical component 
of the electric current density in the photosphere, which is an important measure of the non- 
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potentiality of the active regions, the hehcity content. It has the potential to illuminate 
the active region eruption process. Electric current densities have been repeatedly obtained 



since the photospheric vector magnetic fields became available in 1980s [e.g. ( 


Deloach et 


al. 19841 Lin & Gaizauskas 1987; 


Hagyard||1988 


Canfield et al 


1992 


; Leka et al.| 19931 de 


La Beaujardiere et al. 


1993 Zhang & Wang 1994 


Wang et al. 


1994; 


van Driel-Gesztelyi et 


al. 1994; |Metcalf et al 


1994; Gary & Demoulin 


1995 


Li et al. 


1997; 


Burnette et al. 2004 


Georgouhs & LaBonte||2004 


Gao et al.||2008 


). All these calculations involve the resolution of 



the 180° ambiguity in the transverse field directions (Metcalf et al. 2006). In this paper, we 



will calculate the absolute vertical electric current density, | J^l, with the equation derived by 



Semel & Skumanich (1998) which does not require the 180° disambiguation of the transverse 
field directions. Therefore, this removes one uncertainty in the current calculation. The \ Jz\ 
is particularly useful for studying the mechanical forces due to currents induced in moving 



material (Cowling 1945). 



In section 2, we give a brief description of the IVM instrument and the data reduction 
principles. Sections 3 and 4 correspond to the steps 1 and in the IVM data reduction. 
We describe the procedure used to generate high-quality Stokes images from IVM raw data 
with reference to the active region NOAA 10001 on 20 June 2002 in section 3, and describe 
the procedure to generate the vector magnetic fields from Stokes images in section 4. The 
absolute vertical current densities, | | , are presented in section 5 for two very different active 
regions, NOAA 10001 of 20 June 2002, and 10030 of 15 July 2002. Section 6 gives a summary 
of the results. A detailed description of the magnetic field inversion code, "Triplet Fitter", 
is given in the Appendix. 



2. IVM Data Reduction 

The IVM has a square field-of-view 4.7 arcmin on a side. The basic components of the 
instrument include a polarization modulator, a tunable Fabry-Perot filter, a beamsplitter. 



pre-filters, and Charged-Couple Device (CCD) cameras (Mickey et al. 1996). Among the 
system components, the Fabry-Perot is the fastest at scanning the spectral profiles. The 
time resolution of the system is limited by the CCD camera readout. For the photospheric 
magnetic field observations, the Fabry-Perot filter is centered at Fel A6302.5, with a wave- 
length coverage of 1.2 A and a spectral resolution of 0.04 A. Currently, the IVM detector is 
a 512 X 512 pixel CCD array. Each pixel corresponds to 0.55" x 0.55" spatial size and the 
readout rate per pixel is 320 kHz. It takes 0.8 second to read out the whole array, in which 
time interval is made to a change among the four polarization states, and a scan of a new 
wavelength position. Therefore, an IVM data set including 30 wavelength positions in four 
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polarization states can be completed in ~ 1.6 minutes. 



A special feature of the IVM is that two practically identical cameras record the image 
simultaneously: the so-called "Data" and "Geometry" cameras. Their purpose is to minimize 
the "seeing" effect on the polarization signals obtained in the 4 polarization states at different 
times. The only significant difference between the two cameras is the filters in front of them: 
a narrow band 3 A filter for the "Data" camera to select a single order from the Fabry- 
Perot centered at the Fel A6302.5, and a broad band 100 A filter for the "Geometry" camera 
to select roughly 30 orders from the Fabry-Perot. The idea is that the "Data" camera records 
the polarization images, and the "Geometry" camera records the (nearly) continuum images. 
Any brightness fluctuations in the "Geometry" camera must be due to "seeing" variations, 
while changes recorded by the "Data" camera are due to polarization signals from the Sun. 

A standard IVM calibration package has been implemented in IDL procedures based on 



the IVM observational model described by LaBonte, Mickey & Leka (1999). In the package 



the correction of the dark current, fiat field, scattered light and polarization cross-talk are 
first executed. The effects of the "seeing" on the images include blurring, displacing and 
distorting the fine solar structures. The best "Geometry" image is chosen within a single data 
set as a reference image. Shifting, de-blurring and de-stretching functions are then applied 
to the whole set of "Geometry" images. These very same corrections are then applied to the 



set of "Data" images (Mickey et al. 1996 LaBonte 2004). The final output consists of the 



calibrated Stokes images, /, Q, U, V. 



3. Step 1: Stokes Parameters Reduction 

In this and section 4, we describe the reduction procedures with reference to the active 
region NOAA 10001. In this section, we describe the Stokes parameters reduction corre- 
sponding to the IVM data reduction step 1. On 20 June 2002 this AR consisted of a single 
spot of area 200 [10~^ of solar hemisphere] surrounded by a plage region and projected near 
the disk center N20W07. It was observed continually from 17:00 to 22:20 UT by the IVM, 
giving a total of 149 vector magnetograms. Figjl] shows ARlOOOl on 20 June 2002 taken in 
different seeing conditions in the continuum by the "Geometry" camera (top row), and the 
Stokes I by the "Data" camera at the blue wing (middle row) and the center of the spectral 
line Fel 6302.5 A. 
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17 18 19 20 21 22 

Fig. 1.— Active region NOAA 10001 on 20 June 2002 (N20W07). From top to bottom, 
images are the photospheric continuum taken with the "Geometry" camera (top row), Stokes 
I at the blue wing (middle row) and the Stokes / at the spectral line center (bottom row) by 
the "Data" camera. From left to right, images represent the best seeing, e = 0.0062 at 20:28 
UT, the intermediate seeing, e = 0.0038 at 21:32 UT and the worst seeing, e = 0.0013 at 
22:18 UT. The circles on left-middle image mark the "Quiet Sun" and "Sunspot Penumbra" 
have the radius, r = 2 pixels equivalent to 2" x 2". All images have the same size 240" x 170". 
The plot with circles on the bottom of the figure shows the selected IVM data sets as function 
of time [hours] based on the seeing criterion with seeing condition e > 0.0038. 
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3.1. Data Selection based on the Seeing Condition 

Although it has been best-corrected by the "Geometry" images, the remaining uncor- 
rected "seeing" is the major source of the noise in the calibrated Stokes images. |Leka &: 



Rangarajan (2001) implemented an algorithm for seeing estimation into the IVM software 
package, which is based on the concept of the Modulation Transfer Function. For each 
"Geometry" image, G, the quantity, e, is calculated from root-mean-square (r.m.s.) 

e(A, t) = r.m.s.{[G{X, t) - G,{\, t)]/G,{\ t)} (1) 
where Gg is the image smoothed with a boxcar average. e(A,t), as a function of wavelength 



A at time t, in fact, measures the image contrast. Leka & Rangarajan (2001) commented 
"When this seeing measure decreases towards 0.005, granulation is no longer visible and the 
seeing is estimated 2 arc sec." Ideally, e(A, t) should be invariant with the polarization state. 
In reality, the "Geometry" camera receives the same beam as the "Data" camera which is 
modulated by the polarization modulator. As a result, the Geometry images include the 
polarization characteristics of the spectral line but diluted by a factor of 1/30 because one 
of 30 Fabry-Perot orders falls in the line Fel 6302.5 A. To judge the seeing condition with 
continuum images, we use the values of e(A,t) at the most-blue wing. 

For the 20 June 2002 data, about 20% of the data satisfy e > 0.005. Since we are mostly 
interested in the sunspot and its penumbra, we relax the seeing threshold to e = 0.0038 (see 
the middle column images in Fig. [T]) which allows 53% of the total 149 data sets to be 
selected for further processing. The time distribution of these data sets, which will be used 
to generate the vector magnetograms, is indicated with circles along the horizontal axis in 
the bottom plot of Figjlj The images from the left to right are IVM images with the best 
seeing (e = 0.0062), the seeing value used as the threshold for the current data selection 
(e = 0.0038), and the worst seeing (e = 0.0013). On average, the seeing is ~ 2", which will 
be the minimum pixel sample size in the following data analysis. 



3.2. Spectral Line Shift 

One of the system imperfections is that the central wavelength varies with spatial posi- 
tion on the detector. The central wavelength varies quadratically with position in the field 



because the Fabry-Perot is used in the coUimated beam. This was discussed by LaBonte, 



Mickey & Leka ( 1999 ) and is implemented in the calibration of the fiat field in the IVM soft- 



ware package. The spurious intensity pattern was "simply the convolution of the line profile 
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and the etalon bandpass." The quadratic component of the vignetting is measured in the 
flat-field images taken in the continuum wavelengths. The flat-field image at each wavelength 
is divided by this wavelength-insensitive vignetting. Each fiat field image is fitted again with 
a quadratic spatial distribution which is the wavelength-sensitive quadratic fitting. After 
dividing out the newly fitted wavelength-sensitive quadratic fitting, the fiat-field images are 
multiplied by the continuum quadratic vignetting. 

A time-dependent drift in the central wavelength is also present, but was not discussed 
before. This kind of shift is demonstrated by the Stokes / profiles in Figure |2] For the 
observations with the "Data Camera", we define SDit,x,y, s, X) as the Stokes polarization 
signals (s) at pixel position {x,y) as a function of time {t) and wavelength (A), where s 
represents the Stokes parameters I,Q,U,V. Although the IVM CCD pixel size is 0.55" x 
0.55", the actual spatial resolution is lower because of the atmospheric disturbance and 
residual effect of the atmospheric and instrumental "seeing". The Stokes signals discussed 
here are the average signals over a cluster of pixels. For the current data set, the actual 
spatial resolution is a circular region with r = 2 pixels equivalent to 2" x 2". In Figure [2| 
the Stokes / is calculated as the r.m.s.[SD(t,x,y, s, \)] over the circular "Quiet Sun" region 
shown in Figure [T} After observing for an hour, the amount of the systematic wavelength 
shift is ~ 0.2 A (5 wavelength positions) and, if uncorrected, the spectral line center could 
move out of the spectral passband in only ~3 hours. The wavelength drift is probably caused 
by the Fabry-Perot's transparent plate, although in detail we do not know why, and we do 
not know for how long the wavelength drift has been a feature of the IVM. In the data 



analyzed here, the etalon was periodically re-adjusted (using the method of Mickey (2004)) 
during the observing sequence so that the center of the spectral line re-aligned with the 
center of the Fabry-Perot. This operation is refiected in Fig. [3] in which the Stokes / signals 
at the central position of the IVM spectral window increase with time (top panel) until, 
at UT 19:00 and 21:20, they jump back as the central wavelength position was manually 
adjusted. 

In coincidence, the center of the spectral line varied within the IVM spectral window 
with time (lower panel). We determine the spectral line centers with the minimum Stokes 
/ intensity over the "Quiet Sun". Sometimes, the central positions are one pixel apart in 
consecutive IVM data sets due to the line center fiuctuation. This occurred more frequently 
when the telluric line was strong in the morning such as before 18 UT in the current data 
(see Fig. m. 



The complication caused by the spectral line shift is that it is a mixture of real solar 
vertical motions due to the Doppler effect and the instrumental IVM spectral shifts. The 
uncorrected wavelength drift of ~ 0.2 A in one hour means that, without correction, the 
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Fig. 2. — Examples of the Stokes / spectral line profiles throughout the day. The profiles 
are the signals, r.m.s.[SD{t,x,y, I, \)] over the circled "Quiet Sun" region in Figjl} The 
wavelength shift with time is not uniform in the course of the day. The dips at the wavelength 
position ~ 28 are due to the telluric O2 line. 
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IVM line center drifts by as much as 9 km s~^. The normal IVM integration is 2 minutes, in 
which time the line center drifts by about 0.3 km s~^. The latter sets an absolute lower limit 
to the velocity resolution in the data. It appears that this kind of wavelength shift does not 
vary with the spatial position. 



3.3. Stokes Parameter Uncertainties 

Using the "good seeing" data set, the time-averaged Stokes signals, SD{x,y, s, X), and 
the standard deviations, a{x,y, s, X), represent the statistical averages and uncertainties of 
the single measurement. In Figure|4| Sd{s, x, y, X) (solid lines) and Sd{s, x, y, A)±cr(s, x, y, X) 
(dotted lines) are plotted from two locations, "Quiet Sun", and "Sunspot Penumbra" (see 
circled locations on the middle panel of the first column in the Figjl]). On average, the Stokes 
Q,U,V had uncertainties a = 0.005 in the "Quiet Sun", and a = 0.02 in the "Sunspot 
Penumbra". The signal-to-noise ratios are about 2 for Q, U and 7 for V in both regions. 

The Stokes uncertainties can be reduced by averaging time series IVM data sets. In 
the bottom panel of Figure [3} the IVM data sets are divided into four groups by vertical 
lines. Each group is defined by both the "good seeing" condition and the times at which 
the Fabry-Perot was manually adjusted: 17:01-17:54 UT (27 data sets), 18:03-18:56 UT (15 
data sets), 19:00-21:20 UT (28 data sets), and 21:30-22:11 UT (8 data sets). Within each 
group. Stokes spectral images are averaged at each pixel and wavelength position after (1) 
spatially registering each image to the first one; (2) spectrally registering the spectral line 
to the central IVM spectral window based on the measurements in the "Quiet Sun" (see 
the bottom plot in the Fig. [3|. The signal-to-noise ratios are increased by ~ v^, where 
N = 27,15,28,8, for each averaged Stokes image, respectively. As a result, the signal-to- 
noise levels are enhanced to 10 maximum for Stokes Q, U, and 37 for Stokes V. The Stokes 
uncertainties are 1-2 xlO~'^ in the quiet Sun regions. The four time-averaged Stokes images 
will be used as input to the "Triplet" inversion code for generating the vector magnetograms. 

The time-series averaged Stokes images are highly recommended. Averaging (1) reduces 
the Stokes polarization uncertainties, and increases the signal-to-noise ratios; (2) smoothes 
out rapid fluctuations between the polarization modulations. [These fluctuations are prob- 
ably real on the Sun, but are impossible to be distinguished by the IVM modulation mode] ; 
(3) eliminate oscillatory effects such as the 5 minute oscillation ( Georgoulis fc LaBonte|2006 ). 
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Fig. 3. — Top panel: the Stokes I signals at the central position of the IVM wavelength 
window as a function of time. The signals are taken over the "Quiet Sun" region shown with 
a circle in FigjlJ "*" symbols represent the IVM data sets selected by the seeing threshold. 
Bottom panel: the wavelength positions of the minimum Stokes I intensity as a function 
of time. The "*" symbols are the same data points as those on the top panel. The solid 
vertical lines separate the IVM data into 4 groups. They define the data into four groups 
within which the wavelength shifts were either very small or to the red wing. At ~ 19:12 
UT and 21:30 UT, the central IVM wavelength positions were adjusted back to align with 
the spectral line center by the observer. 
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Fig. 4. — The average Stokes signals and the uncertainties measured at two locations: 
"Quiet Sun" (top four panels); and "Sunspot Penumbra" (bottom four panels) as func- 
tion of wavelength. Solid lines: the time-averaging Stokes signals, SD{x,y, s, X), calcu- 
lated by averaging IVM "good seeing" data sets at each wavelength position; Dotted lines: 
Sjj^x, y, s, A) ± a{x, y, s, A), where a is the standard deviation of all "good seeing" data sets. 
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4. Step 2: Photospheric Vector Magnetic Fields 

This section corresponds to tlie IVM data reduction step 2. Tlie vector magnetic fields in 
NOAA 10001 are inferred from the time-averaged Stokes parameters using the "Triplet" code. 
The latter was developed and implemented in the IVM software package by the late Barry 
J. LaBonte. Because the code has not been described in any published literature, a detailed 
description is given in Appendix A of this paper. In this section, we will briefly describe the 
magnetic field inversion methods in the IVM package. The final vector magnetograms will 
be presented in the second part. The third and fourth sections will demonstrate the effects 
on the magnetograms due to the spectral line shifts with time, and the wrong Stokes V-sign 
set by the IVM hardware. 



4.1. Magnetic Inversion Codes With IVM 



The general problem is to solve the radiative transfer equations for polarized light in the 
solar atmosphere. Many pioneering works have been put forward to lay out the foundation 



for the solutions ( 


Unno||1956; 


Auer, Heasley & House 


1977 


Landolfi & Landi degl'Innocenti 


1982 


Skumanich & Lites 


198' 


7 Jefferies, Lites & Skumanich 1989 


Rees, Murphy & Durrant 



1989). The IVM data reduction package adopts ultimately three methods: (1) the wavelet 
method; (2) the derivative method; and (3) the "Triplet Fitter". Here, we briefly describe 
these three inversion methods. [For the original, but somewhat out-of-date description of 
the magnetic field inversion methods in the IVM package, readers are referred the web site: 



http:/ /www.solar. ifa.hawaii.edu/Reference/IVM/IVM_data_red. html 



The wavelet method, developed by the third author, uses the wavelet transform on 
each profile to locate position and amplitude of the Stokes components. The Paul wavelet 
is used since its real component is similar to Q,U profiles and its imaginary part matches 
the V profile. The line center position is computed from the transform of the / profile, 
then the Q, U, V amplitudes are obtained from their transforms at preselected scales. The 
magnetic field values are obtained by multiplying the polarization parameters by a "magic 
number". Thus, the method suffers from both arbitrariness and likely saturation at large 
field strengths. The magneto-optical effects are not corrected, neither are the magnetic filling 
factors solved. However, this method is quick and, therefore, has been used, for example, 
to reduce a large amount of IVM data for the statistical study of over 1000 fiare-productive 



active regions (Leka & Barnes 2007) . 



The "derivative method" is based on the solution of radiation transfer equations in the 
"weak field" limit, in which the Zeeman splitting is a fraction of the Doppler width (Jefferies, 



Lites &: Skumanich||1989 ) (so-called JLS method). Under the weak field condition, Stokes V 
is proportional to the spectral line slope, dl/dX, and to the line-of-sight field The Stokes 
Q and U are proportional to the derivative of the spectral line slope, d'^I/dX'^, and to the 
square of the transverse field. With a Doppler width for the photospheric line of 40 mA, the 
maximum line-of-sight field strength measured with Fel 6302.5 (Lande factor gj = 2.5) 
is about ~850 G, where Zeeman splitting 4.67 x 10~^^gjX^B\\ A. The method saturates in 
strong field regions and, therefore, is unsuited to magnetic field measurements in sunspots. 
In order to measure strong magnetic fields, an alternative method is needed. 

The "Triplet" code was, therefore, developed as an independent magnetic field diagnosis 



in the IVM package. It originates with the Unno (1956) solutions applicable to homogeneous 



magnetic fields and to absorption coefficients invariant with respect to optical depth. The 



Unno (1956) solutions are incorporated into a non-linear least square method by Auer, 



Heasley & House (1977) to infer the magnetic field from the observed Stokes profiles. The 



"Triplet" code is based on fitting to a more accurate radiative transfer model of the line 



profiles developed by Landolfi & Landi degl'Innocenti (1982) who include the magneto 



optical effect to the non-linear least square method. Another special feature of the "Triplet" 
code is that it treats strong and weak magnetic field regions in two separate steps. In the 
first phase, the photospheric thermodynamic parameters are determined from the continuum 
measurements in the Stokes / spectra. The magnetic parameters are calculated by fitting the 



observed Stokes profiles with equation (A23), which neglects the effects of Faraday rotation. 



In the second phase, the magnetic fields are re-determined in the sunspot regions only, which 
are identified by the relative low brightness in the continuum. The code uses non-linear least 
square fitting of the observed Stokes profiles, including Faraday rotation effects, and initial 
magnetic filling factor < 1. The initial values of the model parameters are those determined 
from the first phase. Most thermodynamic parameters remain unchanged, but some of them 
are allowed to vary. The basis of two-step treatments is demonstrated in the magnetic 



field observations using two magnetically sensitive infrared lines (Lin 1995). By fitting the 



observed Stokes V profile of a normal Zeeman triplet, they found that the photospheric 
magnetic fields consists of two distinct components: the weak field (typically 500 G) is 
found in the intranetwork magnetic elements; the strong field (typically 1400 G) is found in 



the sunspots. Ronan, Mickey & Orral (1987) compared the derivative method and the least 



squares fitting (Skumanich & Lites 1987) of the Unno-Rachkovsky profiles. They concluded 



that the latter is the superior technique for deriving the magnetic field parameters. In this 
paper, we infer the magnetic fields using the "Triplet" code. Before the current paper, the 



code was used to infer the magnetic fields from IVM data by other authors such as Georgoulis 



& LaBonte (2004); Regnier & Canfield (2006); Georgoulis & LaBonte (2006, 2007). 
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4.2. Vector Magnetograms and Uncertainties 

The vector magnetograms are derived from time-averaging Stokes images over 4 time 
intervals: 17:01-17:56 UT, 18:03-18:56 UT, 19:11-21:11 UT, and 21:30-22:11 UT (Figure 
|5|. Before running the "Triplet" code, Stokes- \^ images are flipped in sign. The issue 
will be described in detail in section 4.3. To clearly show the transverse field and the 
magnetic azimuth, the magnetic fields are plotted in size 50" x 50" around the major sunspot: 
the longitudinal fields are the background images; the transverse fields are presented by 
short bars, white bars are against the negative polarities, and black bars are against the 
positive polarities; magnetic azimuths are indicated by the bar directions which remain 180° 
ambiguous. 

The uncertainties of both longitudinal (i?|j) and transverse {B±) fields are estimated 
by calculating the r.m.s.(i?||), r.m.s.{Bj_) over pixels having the Stokes Q,U,V signals less 
than 0.001 representing the quiet Sun. The uncertainties differ slightly among the four final 
vector magnetograms. The errors are 10 G for the longitudinal magnetic component, 40 G 
for the transverse magnetic component, and 0.4 km s~^ for the Doppler velocity. Note that 
the velocity uncertainty is very similar to the spectral line drift within 2 minutes during the 
IVM standard integration time. 

The uncertainties for the magnetic azimuths, 0, are estimated as the r.m.s azimuth 
differences between the time-neighboring magnetograms. The differences are functions of 
the transverse field strengths. The smaller are the transverse field strengths included, the 
larger are the magnetic azimuth uncertainties. When the r.m.s azimuth differences are 
calculated over pixels having 40 < B± < 500 G (representing the plage and sunspot edges), 
the average azimuth uncertainty is ~ 18.5°. When the azimuth uncertainty is estimated over 
the sunspot penumbra (500 < B± < 1500 G), it is on average ~ 8.9°. Figure [6] shows the 
transverse fields of the four magnetograms overlapping one another in the sunspot area. 

4.3. Effects of Spectral Line Shifts 

To quantitatively investigate effects of spectral line shift with time, the four magne- 
tograms described in section 4.2 were used to define the "correct" reference magnetograms. 
The corresponding time-averaged Stokes images described in section 3.3 were used as images 
of "standard Stokes spectra". All the calculations are conducted around the major sunspot 
area size 50"x50", as shown in figure |5} 

We firstly examine differences between the "correct" magnetograms and those derived 
from time- averaging Stokes images without performing the registration of the spectral line 
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Fig. 5. — NOAA 10001 vector magnetograms derived from time-averaged Stokes images 
over four time periods: 17:01-17:56 UT, 18:03-18:56 UT, 19:11-21:11 UT, and 21:30-22:11 
UT. The longitudinal fields are the background images, which strengths are shown in the 
brightness scale [-2500,2500] G. The transverse fields are represented by short bars: white 
bars are against the negative longitudinal polarity, and the white bars against the positive 
polarities. The B± is represented by bar lengths. The distance between neighboring pixels 
represents Bl. = 1500 G. Pixels having Bj_ < 40 G are not plotted. The magnetic azimuths 
are represented by bar orientations. The image sizes are 50" x 50". 
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Fig. 6. — The transverse fields {B± > 500 G) are presented by short bars over-plotting 
the sunspot continuum image. of four magnetograms are distinguished by bar colors: 
17:01-17:56 UT (white); 18:03-18:56 UT (red); 19:11-21:11 UT (yellow) and 21:30-22:11 UT 
(green). When the bar length is equal to the distance between neighboring pixels, B± = 1500 
G. The bar orientations represent the magnetic azimuths, r.m.s. azimuth differences is in 
average ~ 9°. The image sizes are 50" x 50". 
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centers. This was the scenario by Georgouhs & LaBonte (2006) who used time- averaging 
Stokes images over time intervals 20 ~ 30 minutes, and Regnier & Canfield (2006) who used 
the time-averaging Stokes images over time interval 15 minutes. The r.m.s.[AB||, AB±, A0, Av] 
between the two are listed in Table [T] under "Experiment 1" . 

For cases of using non-time averaging Stokes images, i.e., single IVM data set, we 
examine differences between the "correct" magnetograms and those derived from wavelength- 



shifted Stokes images. This was the scenario by e.g. Leka & Barnes (2003a) and Leka & 



Barnes (2007) who used Stokes images over time intervals of 2 or 4 minutes. The wavelength 
shifts are simulated by removing a number of pixels in the red wing from the "standard 
Stokes spectra" images (note that these Stokes images have been time-averaged, therefore, 
the situation is yet different from the real IVM single data set in term of signal-to-noise). 
The removal is equivalent to missing the part of the spectral red wing, but reader should 
note that the number of spectral samples is reduced. The comparisons are listed in the Table 
[l] under "Experiment 2" when the spectral line was shifted by 5 wavelength positions (~ 0.2 
A) to the red wing. This is normally the maximum wavelength shift before the observer 
intervenes to adjust the Fabry-Perot. 

The experiments show that the temporal spectral line shifts do not generally affect the 
magnetic azimuth, neither do they much affect the velocity measurements, but they do affect 
the measurements of the magnetic field strengths. On the other hand, significant departures 
in the magnetic field strengths, magnetic azimuth, and the velocity field from the "correct" 
magnetic field occur when ~ 1/3 line profile has drifted out of the IVM spectral window. At 
this stage, the spectral line center is shifted by ~ 0.44 A to the red wing from the center, 
i.e., leaving IVM running for ~ 2 hours without observer's intervention. On average, the 



r.m.s.ABu ~ 1150, r.m.s.AB± ~ 1024 G, r.m.s.Acp 



14°, and r.m.s.Av 



18 km s 



-1 



4.4. Effects of the V^-Sign Error 



The sign of the Stokes V is set by the polarization modulator in the IVM, but we find 
that the instrumental definition is unfortunately opposite to the conventional definition for 



data taken after January 1999. The standard representation of polarized light (Shurcliff 



1962) implies that Stokes V is negative in the red wing of an absorption line for a magnetic 



field pointing towards the observer ( Rees||1987 Rees, Murphy fc Durrant||1989 ). This follows 



the analysis of polarized radiative transfer (Unno 1956) and is consistent with laboratory 



observations of the Zeeman effect. From the MDI magnetogram (Scherrer et al. 1995), the 
NOAAlOOOl major sunspot polarity is positive. Instead of being negative. Stokes V was 
found positive in the red wing (see the lower-right plot, "Stokes V", in the bottom panel in 
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Table 1. Comparison between the "Correct" and "Incorrect" Magnetic Fields 



Time [UT] ABy [G] ABj_ [G] Av [km s'^] 



Experiment 1 


17:01-17:56 


70.5 


96.5 


5.9° 


0.4 


18:03-18:56 


87.0 


348.4 


13.4° 


0.2 


19:11-21:11 


48.4 


96.4 


6.2° 


0.2 


21:30-22:11 


25.4 


84.9 


9.4° 


0.1 


Experiment 2 


17:01-17:56 


104.9 


84.7 


4.6° 


0.5 


18:03-18:56 


98.3 


125.4 


7.2° 


0.2 


19:11-21:11 


36.4 


76.7 


6.6° 


0.2 


21:30-22:11 


23.2 


62.5 


6.9° 


0.1 


Experiment 3 


17:01-17:56 


1444.9 


337.8 


30.4° 


0.1 


18:03-18:56 


1511.1 


213.5 


26.8° 


0.0 


19:11-21:11 


1507.3 


269.9 


30.3° 


0.0 


21:30-22:11 


1498.8 


169.9 


31.4° 


0.0 



Note. — The numbers are the r.m.s.difFerences between the 
"correct" magnetic fields described in section 4.2 and the "in- 
correct" magnetic fields. The calculations arc carried over pixels 
having the By > 10 and B± > 500 G. Experiment 1: "incor- 
rect" fields arc inferred from time-averaged Stokes images with- 
out spectral line registration; Experiment 2: "incorrect" fields 
are inferred from Stokes spectral images with wavelength shift- 
ing 0.2 A (5 wavelength positions); Experiment 3: "incorrect" 
fields axe inferred from using the wrong Stokes V-signs. 
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Figure |4]). 

If uncorrected in the analysis, the reversed sign of the IVM-^ leads to an over-correction 
of the magneto-optical "Faraday" effect. This is because Faraday rotation occurs when the 
left- and right-circular polarization propagate at different speeds, causing the plane of the 
linear polarization to rotate. Such rotation in the linear polarization is proportional to the 
intensity of the line-of-sight magnetic field component, and is correlated with the sign of 
the component determined by the Stokes V. In a previous version of the IVM software (for 
instance, "Triplet" code), the sign correction was done on the longitudinal magnetic field B\\ 
after the Stokes inversion, so that the sign of B\\ as determined by IVM was consistent with 
results from other instruments. However, this is not the correct approach to the problem 
because it causes inaccurate transverse fields. 

The effects of the reversed Stokes V^-sign is examined in the same way as for the effects 
of the spectral line shifts. The r.m.s. differences are calculated between the "correct" magne- 
tograms and the ones inferred from original IVM Stokes l^-sign. The results are summarized 
in Table [T] under "Experiment 3". The reversal of the V^-sign relative to the conventional def- 
inition not only introduces abundant spurious field azimuths, but also changes the strengths 
of both longitudinal and transverse components of the magnetic field. The Stokes \^-sign 
must be flipped before applying the magnetic field inversion code to the Stokes spectral 
images. 

5. Vertical Electric Current Densities 

The goal of this work is to determine high-quality vector magnetic fields in the photo- 
sphere. In one application, we use the vector magnetic fields to calculate the vertical electric 
current densities for two active regions: NOAA 10001, observed 20 June 2002 UT, was a 
unipolar sunspot region with little change in appearance through the day; and NOAA 10030, 
observed 15 July 2002 UT, was a large quadrupolar sunspot complex showing constant fiux 
convergence, emergence and cancellation during the course of the day. NOAA 10030 pro- 
duced an X3 fiare accompanied by a white light fiare and coronal mass ejections and is the 



subject of several research papers ( 


Liu et. al. 


2003 


Gary and Moore 


2004 


2004 


Li et al. 


2005 


Harra et al. 


2005 


Tian et al. 


2008 


)• 



Georgoulis et al. 
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5.1. Absolute Vertical Electric Current Density 



In SI units, Ampere's law reads //qJ = V x B, where J is the current density and yUo is 
the permittivity of the vacuum. The vertical component of J satisfies 



l^^Jz = (V X B), 



dx 



dy ' 



(2) 



Since /xq = 47r x 10^^ [T m A^^], we can write -iTiJz = (V x B)^, where the magnetic field is 
measured in G and the current density is measured in mA m~^, and the distance is measured 
in m. Determining Jz{x,y) on the photosphere requires knowledge of the direction of the 



transverse field (B^^By), i.e., a solution of the 180° ambiguity problem. However, Semel & 



Skumanich ( 1998 ) showed it is possible to calculate the absolute value of without solving 
the 180° ambiguity problem. They derived the following expression for J'^{x,y): 



(AnJ^y = |(V X B),|^ = B'^g'^ + B^gi - 2{B,By)g,gy, 



where 



Qx = 



9y 



1 

Bi 
1 



^^d{BxBy) 



dy 



^ djBxBy) 
" dx 



2^^^y dy 

2^^^' d^ 



2^^ dx 
2 ^ dy 



(3) 

(4) 
(5) 



From the equations (A26) and (A27), the observable quantities from the IVM are the 
vertical field component, B^ (or the transverse field strength B_i_; and the magnetic 
azimuth, either or + 180°. The two perpendicular horizontal components are written 
Bx = B± COS0 and By = B± sin0, 180° ambiguity in the is equivalent to the observable 
quantity B^By = \B']_ sin 20, and B^ = + By. Neither B^ nor By varies with the and 
+ 180°. In this sense, the IVM provides all necessary quantities for calculating the absolute 
vertical electric current density, | | , which can be determined without solving the ambiguity 
problem. Another condition is that the active region needs to be near the disk center so 
that the longitudinal fields approximate to the vertical fields, Bz, and the transverse fields 
approximate to the horizontal fields, Bj^. 



5.2. Absolute Vertical Current Density in NOAA 10001 



Because the active region NOAA 10001 was near the disk center (N20W07) on 20 June 
2002, the transverse fields are close to the horizontal magnetic fields and the longitudinal 
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fields to the vertical magnetic fields. The absolute vertical current density maps are cal- 
culated from four vector magnetograms (see Fig. [T]). The | J^l noise level is estimated in a 
similar way to the uncertainty of the magnetic azimuth for NOAA 10001. The pixels included 
in the estimation are those having transverse field strengths between 500-1500 G. r.m.s. of 
\Jz\ differences are calculated between neighboring \Jz\ in time sequences, i.e.: |J2|(17:01- 
17:56 UT)-|J2|(18:03-18:56 UT),... 17^1(19:11-21:11 UT)-| 7^1(21:30-22:11 UT), respectively. 
The average noise level is 6.7 mA m~^ for NOAA 10001. 

The strongest currents occur in localized regions at the umbral-penumbral boundary 
where the longitudinal field is about 2000 G, and these currents vary significantly with time. 
More persistent but weaker currents are found in the west and the south-west quadrant of 
the penumbra with magnitudes ~ 15 and ~ 30 mA m~^. These \Jz\ values are well above the 
noise level. Therefore, we believe these signals represent real electric currents on the Sun. 
On the other hand, bright rings surround the spots outside the penumbrae have \Jz\ ~ 10 
mA m~^. They are artifacts of the two phases used in the "Triplet" code, which are amplified 
by the derivatives approximating to finite differences in the \Jz\ calculation. 



5.3. Absolute Vertical Current Density in NOAA 10030 

Figure |8] shows active region NOAA 10030 on 2002 July 15 UT. The spots are marked 
with "PI", "P2", "P3", and "F", for preceding and following sunspots, respectively, in the 
upper-left panel over the MDI continuum. During the course of the day, the major motions 
of the region showed the positive magnetic polarity converging in the directions shown by the 
black arrows in the upper-right panel of the figure. These motions were in close coincidence 
with the emerging magnetic flux, and the flux cancellation around the spot "P3" which also 
rotated counterclockwise. The flare started in the plage area between the filaments Fl and 
F2, and the filament F2 erupted upon flare impulsive onset. 

Vector magnetograms for NOAA 10030 were produced using the same procedure as 
described for NOAA 10001. The IVM raw data were first processed with the IVM standard 
software package. The "seeing" criterion, e > 0.0025, is used for the AR10030, which selected 
~ 40% of the total observed data for further analysis. The Stokes Q,U,V uncertainties are 
2 X 10~^ for a single "good seeing" data set in the quiet Sun. The active region evolved 
rapidly with time. However, time-averaging Stokes images over the scale of an hour is still 
a good choice for increasing the signal-to-noise ratio while maintaining the recognizable fine 
structures at the same time. For all "good seeing" data, the IVM data are divided into three 
groups: 17:30-18:36 UT (29 data sets); 18:36-19:30 UT (18 data sets); and 20:30-22:00 (18 
data sets). Of the three data sets, the first two were taken before a Class X3 flare at 20:03 
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Fig. 7. — The absolute electric current densities are displayed in the background. The super- 
imposed contours represent the longitudinal magnetic field = ±150, ±600, ±1200, ±2400 
G. The "brightness scale" bar in the bottom represents the \Jz\ ranging from to 30 [mA 
m^^]. The image sizes are 50" x 50". 
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Fig. 8. — IVM continuum image image (top) and the IVM longitudinal magnetic field Con- 
tinuum images (bottom) for the NOAA 10030. Both images are taken at 17:35 UT on 
15 July 2002 in the best seeing condition (1"). In the top panel, PI, P2, P3 and F were 
three proceeding sunspots, and the following sunspot, respectively. Two boxes have the size 
50"x50" and 50" x 180", respectively. In the lower panel, Fl and F2 represent two filament 
locations residing above the magnetic field neutral lines. Black arrows indicate the magnetic 
field converging motion in the region. Image sizes are 240" x 150". 
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UT and one of them after the flare. Data taken between 20:00-20:30 UT were not processed 
for magnetic fields because of line distortion by the flare. Nine data sets taken after 22:00 UT 
also satisfied the seeing criterion, but the Stokes spectral lines were badly distorted making 
it difficult for the magnetic field inversion. These data were also abandoned. Three sets of 
time-averaged Stokes images were generated by registering Stokes images in both spatial and 
spectral dimensions within each group. The Stokes uncertainties are 4 — 5 x 10~^. 

Before the magnetic field inversion procedure, the Stokes V signs were multiplied by 
"-" sign. Locating pixels where the Stokes signals < 5 x 10~^, the average uncertainties 
of the magnetic field are 14 G for the longitudinal component, and 44 G for the transverse 
component. The azimuth uncertainty was estimated in the same way as that used for NOAA 
10001, it is 7.5°. Fig|9] shows the three final vector magnetograms from 3 time intervals. 

NOAA 10030 was also near the disk center on 15 July 2002 (N19W01). The transverse 
fields approximate to the horizontal fields and the longitudinal fields approximate to the 
vertical fields. Three maps of the absolute value of the vertical current density, |J^|, are 



shown in Figure 10, \Jz\ uncertainties are estimated in the same way as that estimated for 
the NOAA 10001. \Jz\ uncertainty is estimated as 7.0 mA m~^ which is very similar to the 
uncertainty of NOAA 10001. As an example, we examine the northern area of the "P3" spot 



penumbra which spot rotated ~ 29° a day (see the white arrows in the figure 10) (Li et al. 



2005). The area had average current densities ~ 32 mA m ^ at 17:30-18:36 UT, ~ 24 mA 
m-2 18:36-19:30 UT, and dropped to ~ 18 mA m'^ at 20:30-22:00 UT. These numbers are 
well above the current noise level, implying that the change was real. 



5.4. Relations between Magnetic Field and Current Density 

To investigate the correlation between the current density and the magnetic field, we 
examine the r.m.s.|J^| over pixels having B]\ ± Ai?||. Based on the uncertainty analyses, 
Ai?|| = 5 G. Each magnetogram corresponds to a curve r.m.s.\Jz\ varying with B\\. Figure 
IT] shows relations for NOAA 10001 (top) and NOAA 10030 (bottom), respectively. Instead 
of curves of each magnetogram, the solid curves are the average r.m.s.\Jz\ among four mag- 
netograms for NOAA 10001 (see Fig. [s]), and three magnetograms for NOAA 10030 (see Fig. 
[9]); the dotted curves are the average r.m.s.\Jz\ ± cr (standard deviation of r.m.s.\Jz\ among 
the magnetograms); the dashed lines are the linear least-square fits to the solid curves for 
r.m.s-lJzl above the uncertainty levels. The calculations are carried out over the same areas 
for both active regions, namely 240" x 150". This area is slightly smaller than the images in 
Fig. [l]for NOAA 10001, but is the same as the images in Fig. [8] for AR 10030. 
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Fig. 9. — NOAA 10030 vector magnetograms inferred from three time-averaged Stokes im- 
ages. Backgrounds show the hne-of-sight magnetic fields which strengths are shown in the 
brightness scale on the right, [-2000,2000] G. The transverse fields are represented by short 
bars: white bars against the negative longitudinal polarity, and the black bars against the 
positive polarities. Plotted pixel interval is 6 IVM pixels (0.55" x 0.55" per pixel) which 
corresponds to i?_L = 1500 G. Pixels having B± < 43 G are not plotted. Image sizes are 
240" X 150". 
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Fig. 10. — The absolute vertical current densities are the background images with the NOAA 
10030. The strengths are indicated by the vertical brightness scale on the right, 0-50 [mA 
m~^]. The contours represent the hne-of-sight magnetic field, [-100,100] G. Arrows indicate 
the same region having the | J2I = 32, 24, 18 mA m~^ at different time intervals from the top 
to the bottom plots. Image sizes are 240" x 150". 
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Longitudinal Field Strength (Bn) [G] 



Fig. 11. — r.m.s.{\Jz\) are plotted as functions of the longitudinal magnetic field {B\i) for 
NOAA 10001 on 20 June 2002 (top) and 10030 on 15 July 2002 (bottom). The r.m.s.{\J,\) 
are calculated over pixels having B\\ ±5 G. The solid curves represent the average r.m.s.{\Jz\) 
among magnetograms of four time intervals 17:01-17: 56 UT, 18:03-18:56 UT, 19:11-21:11 
and 21:30-22:11 UT for NOAA 10030; and magnetograms of three time intervals 17:30-18:36 
UT, 18:36-19:30 UT, and 20:30-22:00 UT for NOAA 10030. The dotted curves are the 
average r.m.s.\Jz \ ± a. The dashed curves are the hnear least-square fits to the solid curves. 
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Both active regions show hnear correlations between | and B\\ in the form r.m.s.| = 
a + bxB\\, where a = 3.6856 and b = 0.0044 for NOAA 10001; a = 3.933 and b = 0.0063 for 
NOAA 10030, respectively. Recall the magnetic field force- free field condition, V x B = aB, 
where a is a scalar varying with space. This leads to the vertical component of the equation 
AttJ^ = aB^, where J? is measured in mA m~^, B, ^ B\\ in the current work measured in 



Gauss, and a is measured in m ^. Allowing the force-free field condition valid, the figure 11 
suggests that the vertical current density has a general relation with B\\ 

\j.\ = jo+^m (6) 

where Jo = a, and a/An = b. Jq 7^ 0, but are less than vertical current uncertainties in both 
active regions. 

We should be cautious with how to interpret the equation |6] from observations, which 
seems to suggest that the photospheric magnetic field is linear- force-free, which is an ex- 
treme case requiring a constant-a everywhere in the field. Meanwhile, it is still controversial 
whether the photospheric mag netic field is force-free ( [Metcalf et al]|1995| |Moon et al.1|2002D . 



The figures [7] and 10 clearly show that the a is not constant because there are many fine 



structures in the vertical current density. To further verify non-constant-a within the active 
region, we isolate two areas in the AR 10030, which are marked with boxes on the top image 
in figure [s] The correlation parameters b = —0.0020 (the left box, 50"x50") and b = 0.0031 
(the right box, 50" x 180"), respectively. These local 6s are very different from the general 
b = 0.0063 derived from bigger region 240" x 150". On the other hand, the local b = 0.0040 
which is computed around the major sunspot of NOAA 10001 as those areas shown in figure 
5] (50"x50"). This is very similar to the general b derived with the bigger area 240" x 150". 



Our observations show that \Jz\ increases with increasing By in the photosphere in 
both active regions in large scales, but NOAA 10030 is far from being modelled by a linear 
force-free field. In large scale, \Jz\ increases 1.4 times faster with increasing B\\ in NOAA 
10030 than in NOAA 10001. The ratio of rates of increase varies between 1 and 2 when the 
standard deviations are included in the linear least-square fitting. Recall that NOAA 10030 
produced an X3 flare with emission in the continuum for 6 minutes, and was followed by 
two CMEs. But NOAA 10001 was a simple sunspot region without flare activity. It will be 
interesting to conduct a statistical study using a large sample of active regions. The rate of 
increase may have the potential to forecast the flare productivity of active regions. 

The significant correlations between the total vertical currents and the total magnetic 



flux were reported in a statistical study of over 1000 active regions by Leka & Barnes (2007). 
However, they did not report any differences of correlation parameters among different re- 
gions. It could be that they used an overly simplistic magnetic field inversion-method - 
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"wavelet" method. We also note that the a^i? as the constant-a in the linear force-free 
field was sought by Leka & Skumanich (1999) for a single active region, but no correlation 
was found between the Jz and B\\. They used the "Derivative" method to infer the vector 
magnetic field It would seem prudent to conduct a new investigation of the relation between 
B\\ and \ Jz\ making full use of the current understanding of the Mees IVM. 



6. Summary 

We present the procedures needed to obtain vector magnetic fields from the Imaging 
Vector Magnetograph at Mees Solar Observatory. The general data reduction can be under- 



taken with the IVM software package (LaBonte, Mickey & Leka 1999) but extra steps are 



needed to deal with two issues with IVM and so to generate high quality Stokes images and 
vector magnetograms: 

(1) . Large, previously undiscussed time-dependent spectral line shifts, probably due to 
instability of the Fabry-Perot inside IVM, should be measured. The effects of the wavelength 
shift affect more the magnetic field strengths than the magnetic azimuth or the Doppler 
velocity. 

(2) . A previously undiscussed sign-error built into the IVM definition of the Stokes V 
parameter must be properly corrected. Failure to take this step leads to an over-correction 
for Faraday rotation and to an over-estimated electric current density. In our work, the 
uncorrected l^-sign error introduces average errors A5|| ~ 1490, Ai?_|_ ~ 248 G, and A0 ~ 
30°. Evidence suggests that the sign error occurred after January 1999 when IVM data 
have image sizes 512 x 512. But a careful Stokes- sign verification using the longitudinal 
observations with other instruments is recommended. When Stokes-\^ sign is found to be 
wrong, the Stokes-^ should be flipped by "-" first, then the magnetic fields are inferred from 
the Stokes images. 

(3) Further employing only data taken under good seeing conditions, and correcting 
for the above two errors in the IVM data, we obtained formal uncertainties on the Stokes 
Q, f/, V parameters for the quiet Sun of la ~ 1 x 10"^ for NOAA 10001 and Icr ~ 5 x 10"*^ 
for NOAA 10030. The resulting uncertainties are 10 G for the longitudinal magnetic field 
component, 40 G for the transverse component, and ~ 9° for the magnetic azimuth in the 
sunspot penumbra. 

Absolute vertical electric current densities, \ Jz\ [mA m~^] are calculated for the simple, 
flare-free spot NOAA 10001 on 20 June 2002 UT and the complex active region NOAA 
10030 on 15 July 2002 UT. The calculation does not require disambiguation of the 180° in 
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the transverse field directions. 

(1) The uncertainty on | J,| is ~7.0 mA m^^ for both NOAA 10001 and NOAA 10030. 

(2) In large scale, the vertical current densities, r.m.s.(| Jz|), increase with increasing 
longitudinal magnetic field strength {B\\ ^ B^) in the form | = a + hB\\. The increasing 
rate is 1 - 2 times larger in NOAA 10030 than in NOAA 10001. 

(3) Locally, the linear correlation parameters between |J^| and B\\ are largely variable 
from place to place in complex active region (NOAA 10030), but are fairly constant in the 
simple region (NOAA 10001). 
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work. JL would like to thank David Jewitt for his critical comments throughout the paper, 
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She is also supported by the NASA grant NNG06GE13G to the University of Hawaii to 
support the operation of the Mess Solar Observatory. The work of AvB is supported by 
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A. Magnetic Inversion Method for IVM Data 

Several methods have been developed to infer the vector magnetic field on the Sun 
from polarization measurements of Zeeman-split spectral lines (see section 1). The "Triplet" 
method is based on a simple model for the transfer of polarized light in the solar atmosphere 



(Unno 1956). The model parameters (e.g., magnetic field strength and direction) are varied 



until the best fit to the observed Stokes profiles is obtained (Auer, Heasley & House 1977). 



Here, we use the formulation of the method given by Landolfi & Landi degl'Innocenti (1982) 
(hereafter LL82). We describe the implementation of the method for IVM, which produces 
images in four Stokes parameters (J, Q, t/, V) and 30 wavelengths in Fe I 6302.5 A. The 
implementation takes advantage of the fact that the thermodynamic parameters of the model 
can be estimated directly from the observed continuum intensity. The "Triplet" code was 
originally developed by the late Barry J. LaBonte. 

When a magnetic field is present, the atomic energy levels with angular momentum 
quantum numbers J > are split into multiple components with magnetic quantum numbers 
m = — J, ■ ■ ■ , + J. LL82 described the formation of a "normal" Zeeman triplet in which J = 
for the lower level (no Zeeman splitting) and J = 1 for the upper level. In the case of Fe I 
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6302.5 A, the upper level has J — and the lower level has J — 1, but this reversal of the 
split and unsplit levels compared to a normal triplet does not affect the emergent Stokes 
profiles. Therefore, the formulae presented by LL82 can also be used for the 6302.5 A hne. 

The line formation model used by LL82 makes several assumptions. First, the emitting 
atoms are assumed to be in Local Thermodynamic Equilibrium (LTE), i.e., the populations 

of the atomic levels are in agreement with the Boltzmann equation and the source function 
of the emitted radiation is given by the Planck function, B\(T), where T is the local tem- 
perature. Also, the solar atmosphere is assumed to have plane-parallel stratification in the 
region of emission, and the temperature T{h) is assumed to decrease with height h such that 
the Planck function B{t) at A = 6302 A is a linear function of continuum optical depth, 
B{t) — Bo + BiT (Milne-Eddington approximation). Here Bq is the value of the Planck 
function at the top of the photosphere (temperature minimum region), and Bi is the gra- 
dient of the Planck function. Other parameters of the model are the magnetic field vector 
B, Dopplcr shift AAq (due to mass fiows on the Sun), Doppler width AAd, line broadening 
parameter P, and certain ratios of line- and continuum opacity (see below). These quanti- 
ties are assumed to be independent of continuum optical depth r. The assumption of LTE 
is reasonable in the photosphere where collisional excitation of the atoms dominates over 
radiative processes. The other assumptions are questionable because observations show that 
the solar photosphere is inhomogeneous and dynamic. 

With the above approximations, the equations of radiative transfer of polarized light in 
a Zccman-split spectral line can be solved analytically (LL82). The Stokes parameters of 
the radiation emerging from a magnetic region on the Sun are given by 

/p(AA) = Bo {I + Po^~\i + mm + mf + p']], (Al) 

gp(AA) = -BoPo^-'[{l + m?VQ + {^ + m){'nvPu-r]uPv)+PQXi (A2) 

C/p(AA) = ^Bof^o^-m + mfvu + {l + r]i){r]Qpv-r)vPQ)+Puxl (A3) 

yp(AA) = -Bi/3oA-^[(l + 7y,)V + PFX], (A4) 

where subscript "p" refers to the polarized component of the emission; AA = A — Aq is the 
wavelength offset relative to the rest wavelength Ao; /3o = pBi/Bo; /x = cos^ describes the 
direction of propagation of the light relative to the radially outward direction on the Sun; 
the t]j^Q uy are ratios of line- and continuum opacity; the Pq,u,v describe magneto-optical 
effects; and the quantities A, r;, p and x 3^re defined by 

A = ii + r^,f[{i + r^,f-rf + p']-x\ (A5) 

ri" = VQ + r]l + ril, (A6) 

= Pq + pI + pI, (A7) 

X = VqPq + VuPu + VvPv- (AS) 
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= Kpr-pfc) cos ^, 
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770-^(0, w - t;o), Vb,r = VoH{a, v -vo±vh), 


Pp 




2?7oF(a, V -Vq), pb,r = 2?7oF(a, f - t;o ± t;/^ 



(A9) 
(AlO) 
(All) 
(A12) 
(A13) 
(AM) 
(A15) 

(A16) 
(A17) 

Here v = AX/AXjj, vq = AXq/ AXd and vh = AXh/AXd are the wavelength offset, Doppler 
shift and Zeeman sphtting in units of the Doppler width; a = r/(47rAz/£)) is the damping 
constant; T is the line broadening parameter in frequency units; Auu = cAAd/Aq is the 
Doppler width in frequency units (c is the speed of light); H{a, v) and F(a, v) are the Voigt 
and Faraday- Voigt functions; and t]q is the ratio of line- to continuum opacity in the absence 
of a magnetic field. The Zeeman splitting is given by AXh = 4.6686 x IO^^^qlXqB, where 
B is the magnetic field strength, Aq = 6302.5 A, and ql = 2.5 is the Lande factor for this 
line. The angles and describe the inclination and azimuth of the (constant) magnetic 
field relative to the line of sight. 

The observed Stokes profiles generally include a contribution from instrumental stray 
light. Also, the magnetic field outside sunspots is structured on subarcsecond scales that are 
not resolved by the IVM instrument. Therefore, the observed profiles contain both polarized 
and unpolarized components: 

/(AA) = //p(AA) + (l-/)/s(AA), (A18) 
Q(AA) = /Qp(AA), t/(AA) = /t/p(AA), r(AA) = /Vp(AA), (A19) 

where subscript "s" refers to the unpolarized component of the emission (straylight and/or 
non-magnetic contribution); / is the filling factor, i.e., the fraction of radiation due to 
the polarized component; and /s(AA) is the intensity profile of the unpolarized component 



(Jefferies, Lites & Skumanich 1989). The latter is assumed to be given by equation (Al) 



without a magnetic field: 

/3o 



/s(AA) = fio. 



1 + 



(A20) 
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where Bq s, rjQ s and AAo,s are the source function, hne-to-continuum opacity ratio, and 
Doppler shift characterizing the unpolarization component, and Vq^s = ^^o,s/ ^^d- The 
polarized and unpolarized hues are assumed to have the same Doppler width A Ad, line 
broadening parameter F, and source function gradient /?o- The filling factor / in equation 



( A19 ) varies from pixel to pixel, and cannot be directly determined from the observed Stokes 



profiles. We treat it as a free parameter of the model. Equations (A18) and (A19) give 
the Stokes profile prior to convolution with the instrumental profile. To compare with 
observations we must convolve this profile with the instrumental profile P„(AA), where index 
n indicates the wavelength setting of the IVM instrument (ra = 1, ■ ■ ■ , 30). The profile P„(AA) 
is mostly determined by the Fabry-Perot filter. Hence, the predicted Stokes vector is 

S„ = y S(AA)P„,(AA) dA\, (A21) 



where S(AA) = [/, Q, f/, V] is the Stokes profile given by equations (A18) and (A19). 



The above equations contain 10 unknown parameters {Bq, (3q, rjo, AX^, AAq, F, B, \E', 
(j) and /). To determine these parameters, the "Triplet" code uses two different methods, 
one for plage regions where the magnetic field strength B and filling factor / cannot be 
independently determined (method I), and another for sunspots where a non-linear least- 
square fitting procedure can be used (method II). Method I is also used to obtain initial 
values of the model parameters inside sunspots, as a starting point for non-linear least- 
square fitting. Therefore, we describe this method first. 

Method I: First, certain estimates of the thermodynamic parameters are obtained for 
all pixels in the IVM image. Specifically, the Doppler width is approximated by AXd ^ 
(-^o/c) \/2kT jm + t>|, where T is the brightness temperature associated with the observed 
continuum intensity; m is the atomic mass of iron; and Vt is a micro-turbulent velocity. The 
latter is approximated by Vt ~ [(T/2000) — 1.3] km/s. The line broadening is assumed to 
be dominated by van der Waals broadening, F = 8.08Cg^^M^^n//, where uh = \/2kT /niH is 
the most probable speed of the hydrogen atoms, nu ~ 2 x 10^^ cm~'^ is the total hydrogen 



density at the height where the continuum is formed, and Cg = 8 x 10 (Mihalas 1978). 
The parameters Bq and rj^ are estimated as 

where Ic is the observed continuum intensity, Jq is the observed line center intensity, and 
/?o = 1-3. Next, to obtain more accurate values of the thermodynamic parameters, a quiet 



area in the IVM field of view is selected, and an expression similar to equation (A20) is fit 



to the observed intensity /quiet (A A) from the quiet region. This yields fitted values of the 
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parameters Bq, /Sq, tjq and AXd for each "quiet" pixel. Then the mean of the estimated values 
of Bq for the "quiet" pixels is divided by the mean of the fitted values to obtain a correction 
factor C{Bo), and similar for the parameters f3o, rjo and AAd- Finally, the correction factors 
are applied to the estimated values of Bq, Pq, rjo and AAd for all pixels, including plage and 
sunspot regions (F is unchanged from its estimated value). This yields the final values of 
the thermodynamic parameters for method I. 

We now describe how the magnetic parameters are determined in method I. In facular 
areas outside sunspots, the Zeeman splitting is less than the Doppler width and the Stokes 
profiles do not contain sufficient information to independently determine both the magnetic 
field strength (B) and the filling factor (/); only their product fB is well-constrained by the 
observations. Therefore, we assume a fixed field strength of the magnetic elements, B ^ 1500 
G. Furthermore, we assume that the Stokes profile can be approximated as follows: 

g(AA) = q . ' ^(AA) = u . , V{AX) = v , (A23) 

Sm Wref Sm Wref COS Wref 



where Qrcf(AA) and V^ef(AA) are reference profiles computed with equations (A2) and (A4) 
for some nonzero inclination angle \E'i.ef and azimuth angle ^ref = 0. The parameters q, u 
and V are defined by 

g = / sin^ cos 20, m = / sin^ sin 20, v = f cos'^ , (A24) 



where ^ and are the actual inclination and azimuth angles on the Sun. Equations (A23) 



are fitted to the observed Stokes profiles using linear least-square fitting, which determines 



the parameters g, u and v. Then equations (A24) are inverted as follows: 



t = \/q^^v? = /sin^^, s=^/q^ + u'^ + = f{2 - sin^ ^), (A25) 
/ = |(t + s), \E' = arccos('i;//), = ^ arctan(-u/g). (A26) 

Method II: This inversion method is used in sunspots. The spot areas are identified 
by the differences of the continuum brightness relative to the quiet Sun. For each pixel 
inside a sunspot, the thermodynamic and magnetic parameters are determined by non- 



linear least-square fitting of equation (A21) to the observed normalized Stokes parameters. 



[Q / 1 / 1 / I]n with n = l,---,30. The fitting uses the Levenberg-Marquart algorithm 



(Press et al. 1994). The observed intensities /„ are not used in the fit. Only certain pa- 
rameters are allowed to vary in this fitting process {Bo, rjo, AAq, B, \E', and /), while 
others are fixed to the values determined with method I AXd and F). The reason is 
that the observed Stokes profiles do not provide strong constraints on the latter set of model 
parameters. The initial values of the variable parameters are those obtained from method 
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I. However, to ensure fitting convergence the initial filling factors and field strengths are 
adjusted such that / < 0.8. 

Certain checks are carried out to make sure that the magnetic field strength B > 0, the 
inchnation angle ^ is in the range [0, tt], and the filling factor is in the range / < 1. For 
pixels with parameter values outside these ranges the non-linear fitting process may have 
failed to converge. The parameters at these "bad pixels" are assigned median values from 
surrounding pixels, and the fitting is run again to obtain improved values. The azimuth 
angle is adjusted such that — 7r/2 < < 7r/2 i.e., no ambiguity solution is performed. 

The final output from the "Triplet" code (for both methods I and II) gives the longitu- 
dinal and transverse fields, which are defined by 

S|| = fB cos ^, B^ = fB sin (A27) 

Near solar disk center, B\\ is essentially the magnetic flux density in an IVM pixel. However, 
B± does not have a clear physical interpretation. To obtain the transverse magnetic field 
B± in the unresolved magnetic elements, one should divide B± by the filling factor. 
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